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Can solar wind viscous drag account for CME deceleration? 
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The forces acting on solar Coronal Mass Ejections 
(CMEs) in the interplanetary medium have been evaluated 
so far in terms of an empirical drag coefficient Cd ~ 1 that 
quantifies the role of the aerodynamic drag experienced by 
a typical CME due to its interaction with the ambient solar 
wind. We use a microphysical prescription for viscosity in 
the turbulent solar wind to obtain an analytical model for 
the drag coefficient Cd- This is the first physical charac- 
terization of the aerodynamic drag experienced by CMEs. 
We use this physically motivated prescription for Cd in a 
simple, ID model for CME propagation to obtain velocity 
profiles and travel times that agree well with observations 
of deceleration experienced by fast CMEs. 



1. Background 

The aerodynamic drag experienced by CMEs as they tra- 
verse the interplanetary medium between the Sun and Earth 
is generally thought to arise due to the coupling of the CMEs 
to the ambient solar wind. CMEs which start out slow 
(with respect to the solar wind speed) near the Sun seem 
to accelerate en route to the Earth, while fast CMEs are 
decelerated, suggesting that the solar wind strongly medi- 
ates CME propagation (Gopalswamy et al 2000, Manoharan 
2006) in the interplanetary medium. This fact has been in- 
voked in several papers that derive a heuristic aerodynamic 
drag coefficient for CMEs (e.g., Byrne et al 2010; Maloney 
& Gallagher 2010; Vrsnak et al 2010; 2012; Cargill 2004). In 
particular, Borgazzi et al (2009) have used two different drag 
prescriptions to investigate CME slowdown using a simple 
ID hydrodynamical model that lends itself to analytical so- 
lutions. 

While there has been a fair amount of progress in this 
direction, a physical understanding of the viscosity mecha- 
nism that leads to the drag on CMEs is still lacking. Many 
authors have focussed on the overall dynamics of the CME, 
using the drag coefficient Cd only as an empirical fitting pa- 
rameter that remains constant throughout the propagation 
(e.g., Byrne et al 2010; Maloney & Gallagher 2010; Lara et 
al 2011). 

2. This work 

Our primary focus here is therefore on computing the 
drag force on the CME using a physical prescription. We 
compute the viscosity of the ambient solar wind using a pre- 
scription relevant to collisionless plasmas. We use this so- 
lar wind viscosity prescription to compute the drag on the 
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expanding CME employing approaches that are commonly 
used in fluid dynamics. We restrict ourselves to a simple 
ID, hydrodynamical model in order to focus on the essen- 
tial physics. We first recapitulate the basics of the drag force 
experienced by CMEs. 

3. Viscous drag on a CME 

We start with the oft-used ID equation of motion for 
a CME that experiences only aerodynamic drag (e.g., 
Borgazzi et al 2009): 



mcME Vcme 



rfVcME 
dR 
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where tticme is the CME mass, Vcme is the CME speed, 
Cd is the all-important dimensionless drag coefficient, AT, is 
the proton number density in the ambient solar wind, m p 
is the proton mass, Acme = ""-Rcme is I ne cross-sectional 
area of the CME and V sw is the solar wind speed. The 
term on the right hand side of Eq (1) represents the drag 
force experienced by the CME. It may be noted that this 
form for the drag force is appropriate only for a solid body 
moving through a medium at high Reynolds numbers (e.g., 
Landau & Lifshitz 1987). Most authors adopt an empirical 
drag coefficient Cd that remains constant with heliocentric 
distance. The value for Cd is chosen to ensure that the 
computed velocity profile using an equation like Eq 1 agrees 
with observations (e.g., Lara et al 2011). 

It is well known, however, that the drag coefficient Cd 
is a function of the Reynolds number of the system under 
consideration. We start with data from Achenbach (1972) 
who gives the widely used characterization of Cd as a func- 
tion of the Reynolds number for high Reynolds number flow 
past a solid sphere. This standard characterization can be 
found in most fluid dynamics texts (e.g., Landau & Lifshitz 
1987). The drag coefficient Co(Re) is a very slowly increas- 
ing function of the Reynolds number Re until Re « 10 5 , 
above which it exhibits a sharp drop; this is the so-called 
drag crisis. Beyond this sharp drop, Cd is an increasing 
function of Re; this spans the supercritical and the trans- 
critical regimes. While these facts are well known, we are not 
aware of any analytical fits to Achenbach's (1972) data, es- 
pecially for high Reynolds numbers. We therefore fit Achen- 
bach's (1972) data for Re > 10 6 to the following functional 
form: 

Cn(Re) = 0.1478- +9.8 x 10" 9 .Re. (2) 

Re 

In our case, the Reynolds number is a function of the 
CME velocity Vcme, the typical macroscopic lengthscale 
Acme and the viscosity coefficient according to the usual 
formula 

„ Vcme -Rcme Vcme Rcme Ni m p 

Re = = , (6) 

v 7] 

where v (cm 2 s _1 ) is the coefficient of kinematic viscosity 
and r\ = vNi m p . Our use of the CME radius Rcme for a 
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typical macroscopic lengthscale implies that (as for a flux 
rope CME) the cross-sectional dimensions of the CME are 
~ -Rome. In what follows, we will develop a physically mo- 
tivated viscosity prescription for the solar wind and use it 
to determine the Reynolds number self-consistently. We will 
use the Reynolds number thus determined to calculate Cd 
and use it to solve the CME equation of motion (Eq 1). 

4. Viscosity prescription for a collisionless 
plasma 

It is well known that the solar wind is a collisionless 
plasma above ~ 10 Rq, where the (Coulomb) collision mean 
free path is longer than macroscopic scale lengths. How- 
ever, a fluid description is a generally well accepted one for 
the solar wind. This implies that the collisionless particles 
are confined via effective collisions with scattering centers of 
some sort, which lends validity to a fluid treatment. Clearly, 
viscosity in the solar wind is not due to interparticle colli- 
sions, as is the case in everyday experience. 

4.1. Hybrid viscosity 

We follow Subramanian, Becker & Kafatos (1996; SBK96 
from now on), who developed a model for collisionless vis- 
cosity. They considered a one-dimensional, plane-parallel 
shear flow, with viscosity being provided by the flux of pro- 
tons that originate from one of the layers and impinge on 
the other. Tangled magnetic fields are embedded in the 
flow, and the average coherence length (i.e., the length for 
which an average magnetic field line is expected to remain 
straight) of the "kinks" in the magnetic fields is taken to 
be A co h- The protons, whose gyroradii are negligible in 
comparison to other length scales, are envisaged to slide 
along the field lines, and change their direction (i.e., get 
scattered) either when they encounter another proton, or a 
kink in the tangled magnetic field. Since the protons are 
collisionless, proton-proton collisions are unlikely, and mo- 
mentum transfer (and consequently viscosity) is dominated 
by proton-magnetic kink encounters. The magnetic kinks 
can also be envisaged as turbulent eddies that act as scat- 
tering centers; the only restriction in this scenario is that the 
scattering centers evolve over a timescale that is slower than 
that of the travel time of a typical proton, so that they ap- 
pear stationary to it. Using this scenario, SBK96 arrive at 
the following simple and physically motivated formula for 
a "hybrid" coefficient of dynamic viscosity rfryb! one that 
is neither due to magnetic field stresses alone, nor due to 
proton-proton (Coulomb) collisions alone: 



Vhyb = Vhyb Ni m p = — — r/s g cm 1 s 1 , (4) 
15 An 

where A denotes the effective mean free path, Ah is the 
(Coulomb) mean free path for proton-proton collisions and 
rjff is the standard proton viscosity due to Coulomb colli- 
sions alone (Spitzer 1962). The viscosity ?7hyb in a collision- 
less plasma is thus suppressed with respect to the standard 
Coulomb value ijg by a factor (2/15)A/Aii. For the sake of 
completness, we reproduce the expressions for Ah and t]b- 

Ah = twin = 1.8 x 10 5 -— 7— - cm 
Ni InA 

j,5/2 

r)g — V&Ni m p iirms An = 2.2 x 1CT 1 '' ' — gem -1 s _1 (5) 

where v rms = (SkTi/mp) 1 ^ 2 is the rms thermal velocity of 
the protons, tu is the mean interval between proton-proton 
(Coulomb) collisions, Ni is the proton number density in 



cm -3 , Ti is the proton temperature in Kelvin and InA is 
the Coulomb logarithm, which is taken to be equal to 20. 
Eqs (4) and (5) comprise our operational definition for hy- 
brid viscosity, which we will apply to the solar wind. It may 
be noted that the hybrid viscosity (Eq 4) can be equivalently 
written as 

2 

?7hyb = — V6 Ni m p «rms A . (6) 

15 

Eq 6 looks similar to the usual fluid dynamical expression 
for turbulent viscosity (e.g., Landau & Lifshitz 1987) 

??turb ~ Ni trip Vturb hurb , (7) 

where v tU ib is the velocity of the turbulent eddies of length- 
scale /turb- Furthermore, if hurb is interpreted as the dissipa- 
tion lengthscale in a turbulent cascade, and if we adopt the 
usual expression for energy dissipation rate in Kolmogorov 
turbulence e ~ v 3 urb /Zturb, the expression for turbulent vis- 
cosity can be rewritten as (Verma 1996) 

Jfturb ~ Ni m p (e/turb) 1/3 ■ (8) 



4.2. Solar wind viscosity 

We now use Eqs (4) and (5) to compute the operative 
viscosity in the ambient solar wind. As mentioned earlier, 
the viscosity would be determined primarily by collisions 
between protons and magnetic kinks, and therefore the ef- 
fective mean free path A « A co h- The coherence length of 
the magnetic field irregularities (A co h) represents the short- 
est lengthscale over which the turbulent magnetic field is 
structured. It is fairly well established that density turbu- 
lence in the solar wind follows the Kolmogorov scaling, with 
an inner (dissipation) scale that is determined by proton cy- 
clotron resonance (Coles & Harmon 1989). We assume that 
the magnetic field irregularities follow the density irregular- 
ities (e.g., Spangler 2002), and are governed by the same 
inner scale. We therefore take A co h to be equal to the inner 
scale of solar wind density turbulence advocated by Coles & 
Harmon (1989): 

A^A coh = 684A^ 1/2 km. (9) 

The solar wind proton density N is assumed to be given 
by the model of LeBlanc et al (1998): 

Ni(R) = 3.3 x 10 5 R~ 2 + 4.1 x 10 6 R~ 4 + 8 x 10 7 R~ 6 cm" 3 , 

(10) 

where R is the heliocentric distance in solar radii. We take 
the proton temperature to be Ti = 10 5 K. 

Thus Eq (4), with A given by Eq (9) and the density given 
by Eq (10) defines the hybrid viscosity prescription for the 
ambient solar wind for our purposes. This enables us to ob- 
tain the viscosity of the ambient solar wind as a function of 
heliocentric distance. 

Eviatar & Wolf (1968) obtain an estimate of v « 
800 km 2 s -1 at the Earth by considering the momentum 
transfer across the magnetopause due to fluctuations in- 
duced by the two-stream cyclotron instability. This is in 
excellent agreement with the value of 788 km 2 s _1 we ob- 
tain for the kinematic viscosity coefficient from the hybrid 
viscosity model. Perez-de-Tejada (2005) has derived a rough 
estimate of v ~ 1000 km 2 s~ for the kinematic viscosity of 
the solar wind near the ionosheath of Venus (heliocentric dis- 
tance 0.72 AU). The hybrid viscosity model yields a value of 
600 km 2 s -1 for the coefficient of dynamic viscosity at 0.72 
AU. 
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5. Results 

We use the prescription for the viscosity given in § 4 to 
determine the Reynolds number (Eq 3), which in turn is 
used to determine Cd (Eq 2). The resulting expression for 
Cd is used in the equation of motion (Eq 1) to solve for the 
CME speed as a function of heliocentric distance. 

A representative result is shown in Figures 1 and 2. In 
this example, the CME radius -Rome is assumed to expand 
as -Rome = KRP , where K is a constant of proportionality, 
which we determine by assuming that the radius of the CME 
is 1 Rq at a heliocentric distance of 2 Rq. The power law 
index p is assumed to be equal to 0.78 (Bothmer & Schwenn 
1998). The solar wind speed V svl is taken to be 375 km/s 
and the velocity of the CME at R = 5Rq is taken to be 
equal to 1048 km/s. This CME initial velocity is represen- 
tative of a CME on May 17 2008 that was well observed by 
the STEREO spacecraft from the Sun to the Earth (Wood 
et al 2009). We chose this CME since it is one of the few fast 
ones that were well observed from the Sun to the Earth by 
the STEREO spacecraft during the minimum and ascend- 
ing phases of cycle 24. The CME mass tocme is taken to be 
5 x 10 14 g. This value for the CME mass is representative 
of CMEs during this part of the solar cycle (see Vourlidas 
et al 2010). 

The solid line in Fig 1 is the solution for the velocity pro- 
file calculated from Eq (1). The plus signs represent data 
for the CME of May 17 2008. In order to get the velocity- 
distance data for this CME, we started with the distance- 
time data for this CME (top panel, Fig 6 of Wood et al 
2009). The velocity-distance data points in Fig 1 are derived 
by numerically differentiating the distance-time data. The 
large undulations in the velocity-distance data shown in Fig 
1, especially for distances > 20_Rq, are due to fluctuations 
in the distance-time measurements, which are accentuated 
in the derived velocity. In particular, the significant dip in 
CME speed after ~ 150 Rq is almost certainly an unphysical 
artefact. It is well known that CMEs achieve an asymototic 
speed by the time they reach the Earth, and often as soon 
as ~ 100 Rq (e.g., Poomvises, Zhang & Olmedo 2010). The 
dip in the CME speed beyond 150 Rq is probably due to the 
difficulty in tracking CME in HI data, which are well known 
to have poor signal to noise ratio at those distances. In 
view of this, the agreement between the solid line, which is 
the solution to Eq (1) and the velocity-distance data for the 
May 17 2008 CME is evidently rather good. The Reynolds 
number for the CME (in units of 10 s ) and the dimensionless 
constant Cd (which is now computed self-consistently, using 
the prescription for solar wind viscosity) are shown in Fig 2. 
The results shown in Fig 2 justify the function we adopt for 
Cd (Eq 2), which is valid only for Reynolds number S> 10 4 . 
We note that Wood et al (2009) obtain the velocity for the 
May 17 2008 not by piecewise numerical differentiation of 
the distance-time data (as we have done), but by assuming 
an empirical model that incorporates an initial acceleration 
phase and a subsequent deceleration phase that is followed 
by a constant velocity phase. 

We have also carried out this exercise for some of the 
fastest Earth-directed halo CMEs observed so far. Our 
results are summarized in Table 1. For each event, 
the parameters supplied to the model (Eq 1) are the 
solar wind speed Vsw, the CME mass hicme and the 
CME initial velocity Vt. The initial velocity of each 
halo CME is determined from the LASCO CME catalog 
(http://cdaw.gsfc.nasa.gov/CME_list/). Since these are 
halo CMEs, their masses are hard to estimate, and we have 
used reasonable guesses, based on typical CME masses dur- 
ing the appropriate phase in the solar cycle (Vourlidas et al 



2010). As before, the proton temperature is assumed to be 
10 5 K and the Leblanc et al (1998) density model (Eq 10) is 
used for determining the solar wind viscosity and the drag 
coefficient Cd- The quantity Vicme represents the speed of 
the relevant interplanetary CME detected in-situ by space- 
craft near the Earth and TTicme is the time elapsed between 
the first detection of the halo CME in the LASCO FOV 
and the detection of the corresponding ICME at the Earth. 
The quantity Vmodel represents the predicted ICME velocity 
at the Earth and TT mo d e i denotes the predicted Sun-Earth 
travel time. It is evident from Table 1 that the model predic- 
tions for the ICME speed at the Earth and the total travel 
time agree quite well with the observations. 

We have used the first CME in Table 1 (20010409) as a 
representative event to ascertain the effects of varying some 
of the parameters on the model predictions. All other quan- 
tities remaing fixed, we find that a 10 % decrease (increase) 
in \4w results in a 7.5 % decrease (increase) in Vmodd and a 
7 % increase (decrease) in TTmodei- Similarly, we find that 
a 50 % decrease (increase) in racME leads to a 13 % de- 
crease (increase) in Vmodel and a 15 % increase (decrease) in 
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Figure 1. The solid line denotes the model prediction 
for the CME velocity (in km/s) as a function of heliocen- 
tric distance (in Rq). The plus signs are derived from 
the distance-time data for the May 17 2008 CME. Fur- 
ther details are mentioned in the text. 
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Figure 2. Heliocentric variation of the CME Reynolds 
number in multiples of 10 8 (RN) and the drag coefficient 
(Cd) for the model of Fig 1. 
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Table 1. A comparison between observations and model results for near-Earth ICME speeds and Sun-Earth travel times for some 
fast Earth-directed CMEs. 

Event Vsw "cme Vj Vicme V'modci TT ICME TT modcl 
(km/s) [g] (km/s) (km/s) (km/s) hours hours 

20010409 450 5 x 10 14 1192 670 696 53A 533 
20031028 450 2 x 10 15 2459 1400 1350 24.8 25.6 
20040726 400 10 15 1366 900 894 38.63 37 



TTmodci- On the other hand, a 100 % increase (decrease) in 
Ti results in a 4 % increase (decrease) in Knodei and a 4.6 
% decrease (increase) in TT m odei- Therefore, the results of 
our model are most (least) sensitive to the solar wind speed 
(proton temperature). 

6. Summary and Conclusions 

6.1. Summary 

To summarize, we have considered the motion of a CME 
under the influence of only a drag force (Eq 1). The drag 
force is proportional to the square of the CME speed (rela- 
tive to the background solar wind); this form is appropriate 
for a solid body moving through a background medium at 
high Reynolds numbers. The proportionality constant in- 
volves the important (dimensionless) drag coefficient Cd, 
which has been treated as an empirical fitting parameter in 
the literature so far. 

It is, however, well known that the drag coefficient Cd 
is in fact a function of the CME Reynolds number Re. We 
use standard data pertaining to the motion of a solid sphere 
at high Reynolds numbers to characterize Cd as a function 
of Re (Eq 2). It is necessary to know the coefficient of 
kinematic viscosity v in order to determine Re (Eq 3). We 
compute the operative viscosity in the (quiescent) collision- 
less solar wind using a prescription that considers protons 
colliding with magnetic scattering centers (Eqs 4 and 5). 
We assume that the proton temperature in the solar wind is 
T; = 10 5 K, the proton number density is given by Eq (10) 
and that the operative mean free path A is the inner scale of 
the spectrum of density irregularities (Eq 9). This gives a 
simple, physically motivated prescription for the drag coeffi- 
cient Cd, which we use in solving a simplified ID equation of 
motion (Eq 1). Using observational estimates for the CME 
starting speed and reasonable guesses for the ambient solar 
wind speed and the CME mass, we find that our results for 
the CME speed profile and the Sun-Earth travel time agree 
quite well with observations (Fig 1 and Table 1). 

6.2. Conclusions 

The agreement between our theoretical predictions and 
observations of CME deceleration is remarkable, especially 
in light of the several simplifying assumptions we have 
adopted. It confirms the essential validity of our physical 
prescription for the dimensionless viscous drag parameter 
Cd, which can be profitably used in semi-analytical treat- 
ments (e.g., Vrsnak et al 2012) and simulations of CME 
propagation. 

We next mention some caveats that accompany our con- 
clusions. Firstly, this is a rather simplified, ID hydrody- 
namic treatment that only addresses the essential physics 
of the CME-solar wind interaction. There is no attempt to 
include Lorentz force driving of CMEs; something that is 
known to be important upto around 30 Rq (e.g., Subrama- 
nian & Vourlidas 2007), at least for CMEs that are only 
moderately fast. Very fast CMEs, such as the ones consid- 
ered in Table 1, probably experience Lorentz force driving 
and acceleration very early in their evolution. Upcoming 
instruments such as the ADITYA-I coronograph (Singh et 



al 2011) could address this issue via high cadence images of 
the inner solar corona. 

Secondly, it is surprising that a drag law of the form of 
Eq (1), which is in fact valid only for solid bodies, works so 
well. CMEs are often thought of as bubbles, which are techni- 
cally defined as bodies that deform in response to tangential 
stresses on their surfaces. The drag force for high Reynolds 
number flow past a bubble is in fact proportional to Vcme, 
and not to Uc ME (Landau & Lifshitz 1987; Merle et al 2005). 
With this form for the drag force, our preliminary results in- 
dicate that a typical fast CME slows down only by about 1 % 
by the time it reaches the Earth. It is therefore worth inves- 
tigating why CMEs seem to behave like solid bodies as far as 
their interaction with the ambient solar wind is concerned. 
In a flux rope CME, the large-scale, ordered magnetic field 
of the flux rope could provide an explanation for this ap- 
parent "solid body" effect. Finally, we note that the drag 
coefficient derived from ideal MHD simulations (e.g., Cargill 
2004) can often significantly exceed unity, in contrast to the 
values obtained in this work. 
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